Diffusion Resonances in Action Space 
for an Atom Optics Kicked Rotor with Decoherence 
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We numerically investigate momentum diffusion rates for the pulse kicked rotor across the quan- 
tum to classical transition as the dynamics are made more macroscopic by increasing the total 
system action. For initial and late time rates we observe an enhanced diffusion peak which shifts 
and scales with changing kick strength, and we also observe distinctive peaks around quantum reso- 
nances. Our investigations take place in the context of a system of ultracold atoms which is coupled 
to its environment via spontaneous emission decoherence, and the effects should be realisable in 
ongoing experiments. 

PACS numbers: 05.45.Mt, 03.65.Yz, 42.50.Lc 



The transition from quantum to classical behaviour in 
nonlinear dynamical systems has been a topic of much in- 
terest in recent years, motivated by the large differences 
that exist between the behaviour of such systems in the 
two regimes. The 5-Kicked Rotor(DKR) is a prime ex- 
ample - in particular, coherence effects in the quantum 
DKR completely suppress classical diffusion [jl], p|. The 
quantum DKR is also very interesting because it has been 
beautifully demonstrated in experiments which probe the 
momentum distribution of a cloud of laser-cooled atoms 
interacting with a pulsed standing wave of near resonant 
light H. These experiments necessarily involve a gener- 
alisation of the DKR to kicks of finite length, and it is 
this system, the Kicked Rotor(KR) which we investigate 
in this paper. 

There have been numerous studies, both theoreti- 
cal (see for example 0, ||, ^, |§]) and experimental 
||, [l0| [ll], [l^, |l3| , of the role of decoherence in the quan- 
tum to classical transition for the Kicked Rotor. As with 
any real quantum system, the atom optics KR couples to 
its environment, resulting in a loss of phase coherence. 
In the case we consider here, this coupling is between 
the atoms and the vacuum electromagnetic field, and re- 
sults in atomic spontaneous emissions and concomitant 
random momentum recoils. 

Most work in the past has focussed on changing the 
level of this decoherence and observing the effects on 
momentum diffusion rates and distributions, looking in 
particular at how increased levels of decoherence "drive" 
the quantum system towards classical behaviour. In this 
paper we focus instead on what happens when we fix 
the level of decoherence in our real quantum system, and 
then make the dynamics more macroscopic by varying 
the total action in the system; that is, by varying the 
effective Planck's constant. In so doing we find dramatic 
structures in the momentum diffusion rates, similar to 
those found recently for the DKR with a continuous po- 
sition measurement |l4|] , only here we consider a specific 
experimental configuration. 

The system we model is a cloud of ultracold Caesium 



atoms (of initial temperature s=s lOpK ) which interact 
with a standing wave of laser light of frequency loi, de- 
tuned far from the frequency wq of the 6S1/2 — > 6P3/2 
atomic transition. The laser is pulsed with period T and 
pulse profile f(t). If the detuning is large, the internal 
atomic dynamics can be eliminated, and the resulting 
single particle Hamiltonian (for just the external degrees 
of freedom) is || 
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where x and p are the atomic position and momen- 
tum operators, respectively, and ki is the wave num- 
ber of the laser light. The effective potential strength, 
Q e ff = ri 2 (s45/^45 + S44/644 + 543/543), accounts for the 
different dipole transitions between hyperfine levels in 
the Caesium atoms (F=4 — * F'=5,4,3), with 8ij the cor- 
responding detunings, and Q/2 the (single-beam) reso- 
nant Rabi frequency. If we assume equal populations 
in all Zeeman sublevels, then S45 = S44 = and 



S43 = y^g . We can rewrite this Hamiltonian in appropri- 
ate dimensionless units as 
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which is the Hamiltonian for the standard kicked rotor 
system. Here, <j> — 2kix, p = 2kiTp/m, t' = t/T, and 
H' = (AkfT 2 /m)H. The classical stochasticity param- 
eter is given by n — Q e j fU)nTT p , where t p is the pulse 
length and u)r — hkf /2m. In our work f(t') is gener- 
ally a square pulse, i.e., f(t') = 1 for < t' < a, where 
a = Tp/T. Note that k — n/a. In these units, we have 
[<fi, p] = ik, with k = &ojrT. Thus the quantum nature 
of the system is reflected by an effective Planck's con- 
stant, k, which scales as we change the total action in 
the system by altering the pulse period T. 

Decoherence occurs in the form of spontaneous emis- 
sion events, which occur when the atoms absorb light 
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from the standing wave M. It is assumed that momen- 
tum distributions in orthogonal directions remain inde- 
pendent, and thus the system remains effectively one di- 
mensional. We characterise the level of this decoherence 
by the probability of spontaneous emission per kick, rj. 
Given the large detuning, i.e., Sl e ff/S <C 1, this process 
may be modelled by the master equation for the density 
operator w of the system M 
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where N{u) is the distribution of recoil momenta pro- 
jected onto the standing wave axis, and [., .]+ denotes an 
anti-commutator. We have also modelled spontaneous 
emission events in which the atoms absorb light from 
counterpropagating beams of oppositely circularly po- 
larised light, which interact with the atoms continuously 
p0| . This leads to emission events which are indepen- 
dent of position, and which may occur at times other 
than during a kick. The results obtained for the two 
types of spontaneous emission noise are very similar, and 
so we only present the results for the first type here. 
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FIG. 1: Momentum diffusion rates D(n) in the quantum 
kicked rotor as a function of the effective Planck's constant, 
R for (a) the first two kicks (n = and n = 1) and (b) the 
third, fourth and fifth kicks (n = 2, n = 3 and n = 4), with 
r\ = 10%, k = 9, and a — 0.005. Classical values are marked 
on the vertical axis. 

We simulate Eq. (0) using the method of quantum tra- 
jectories, as described in Ref. B. We choose initial 



momentum eigenstates from a Gaussian distribution of 
width <7 p /k = a p /2hki = 4 (which reflects the initial tem- 
perature distribution) and we take an incoherent average 
over the final momentum distributions. We choose a to 
be small so that the effects of KAM boundaries are not 
important Typically we use 1000 trajectories, and 
we calculate statistical errors (shown for some points in 
the figures) based on dividing these trajectories equally 
into 10 groups and computing errors in the means. We 
are primarily interested in the momentum diffusion rate, 
which is defined as the change in the kinetic energy from 
one kick to the next, 2D(n) =< f>n+i > ~ < Pn >> where 
we denote po = p(t' = 0), p\ = p(t' = 1) etc. 

Our simulation results for the momentum diffusion 
rates across the first five kicks are shown in Fig. 0. Aside 
from some noise in the simulations for larger values of k, 
we see that for the first two kicks the diffusion rates are 
essentially constant with respect to S. This is the quan- 
tum version of the quasilinear behaviour known to exist 
classically in the kicked rotor system. In fact one can 
show (for the DKR) that for a uniform initial distribu- 
tion of positions, 2D(0) = n 2 /2 (shown as a solid line in 
Fig. 0(a)). Similarly, if we also assume an initial Gaus- 
sian momentum distribution of standard deviation a p , it 
can be shown that 
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where J n is an ordinary Bessel function of order n, K q = 
2Ksin(£/2)/£ and K 2q = 2Ks'm(k)/k. For sufficiently 
large tr p , this reduces to the same result as for the first 
kick, which can be seen for our system in Fig. 0(a). 

After the second kick, the system settles down into 
its initial quantum diffusion period, where for a small 
time the system exhibits classical-like diffusion, with a 
relatively constant momentum diffusion rate, D q . As can 
be seen from Fig. 0(b) , the k dependence of this rate is 
quite remarkable. We observe an enhanced diffusion peak 
(or resonance) around k — 3 which shifts to the right 
and increases in magnitude as we increase k (See Fig. 
||). There is also a peak in the diffusion rates near the 
quantum resonance at k = 2ir [ p^[ . (This structure also 
occurs at larger multiples of 27r) . The quantum diffusion 
rates that we observe in this regime agree well with those 
predicted by Shepelyansky under the conditions % > 1 
and k^$> k lltj, i.e., 



HK q ) - Jl{K q ) + J 2 {K q ) + 4(K q \ 
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This can be seen in Fig. |2|, where we plot the average 
of the curves for D(2 — 5). There is surprisingly good 
agreement for the quantum resonance peak and for the 
position of the enhanced diffusion peak, especially consid- 
ering that the condition k>S does not hold for our large 
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k values. The discrepancy in the height of the peak is cre- 
ated mainly by our choice to average over the diffusion 
rates from 4 different kicks. It is a qualitative decision 
as to when the system has really settled into the initial 
quantum diffusion regime, but averaging over D(2 — 5) as 
we do gives us an objective estimate of the corresponding 
diffusion rate. However, in some cases (particularly near 
the maximum of the enhanced diffusion peak), the dif- 
fusion rate will already have begun to decrease towards 
the late time diffusion regime before the sixth kick. 

The initial quantum diffusion period lasts for a small 
number of kicks, after which the diffusion rate begins to 
decrease. In the absence of noise, the system settles into 
a localised state |3j, and D(n) — > as n — > oo. However, 
the onset of dynamical localisation is a coherence effect, 
and in the presence of decoherence the system settles into 
a late time diffusion regime where D(n) — > Doo ^ as 

71 — » OO. 

Our simulation results for these late time diffusion 
rates are shown for varying levels of decoherence, 77, in 
Fig. |[ and for varying k in Fig. |i[ Again we observe 
an enhanced diffusion peak (or resonance) which shifts 
and scales with increasing k, as well as a much more nar- 
row peak near the quantum resonance at % = 2n. The 
most notable feature here is that for appropriate values 
of k and 77, the momentum diffusion rates near the top 
of the enhanced diffusion peak are actually larger than 
the corresponding classical values (marked on the graph) , 
even if we account for added momentum diffusion due to 
spontaneous emission in the classical system. 

The late time diffusion rates may be approximated by 
the formula |j| |) 

00 

A» = 5>(1 - »?) n A>(n), (6) 
n=0 

where Dq (n) is the diffusion rate at the nth kick for a KR 
without decoherence. Essentially, we assume here that 
dynamical correlations over particular time intervals are 
suppressed by a factor expressing the probability that 
a spontaneous emission occurs within that time interval. 
The correlations taken over a set number of kicks give rise 
to the diffusion rates seen in the KR without decoherence 
after that number of kicks. Hence, we take a weighted 
average over the diffusion rates as the KR settles down, 
where the weighting for Do(n) gives the probability that 
the first spontaneous emission event occurs during kick 
number n + I, i.e., 7/(1 — 77)™. The early time diffusion 
rates are thus "locked in" by the loss of phase coherence. 

In Fig. ^ the simulation results for are plotted as 
points, and calculations of the right hand side of Eq.(||) 
based on calculations of Do(n) from simulations of the 
KR without decoherence are shown as solid lines. There 
is very good agreement between the two sets of values, 
especially for the enhanced diffusion peak. Note that the 
values given by the solid lines contain statistical errors 
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FIG. 2: Initial quantum diffusion rates in the kicked rotor 
for varying values of k and k, with r\ = 10% and a — 0.005. 
The points show simulation results (an average over D(2 — 5)), 
while the lines show Shepelyansky's formula, Eq.(^). Classical 
values are marked as points for k = 0. Note the use of a linear 
scale for k. 
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FIG. 3: Late time momentum diffusion rates in the quantum 
kicked rotor for varying values of rj and k, with k = 10 and 
a = 0.005. The insert shows the peak near the quantum 
resonance at k = 2tt, and classical values are marked for k = 0. 
Note the use of a linear scale for k. 

from the simulations of Dq(ji) which are comparable in 
magnitude to those displayed in the figure. The level 
of agreement indicates that the model associated with 
Eq.(^J) works very well for the late time diffusion rates. 

Late time diffusion rates greater than the correspond- 
ing classical rates occur because quantum correlations 
cause the initial quantum diffusion rates to be higher 
than the corresponding classical rates in the appropriate 
cases, and these higher rates become locked in by the loss 
of phase coherence. These decoherence effects are thus 
much more important than the more direct increase in 
momentum diffusion due to the recoil in a spontaneous 
emission process which may be accounted for classically 
as well as quantum mechanically. 

It is possible to find an analytical expression for 
which agrees well with the simulation results over a 
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FIG. 4: Graph showing simulated late time momentum diffu- 
sion rates (points) and results from Eq.(|^) (solid lines), with 
r\ = 10% and a — 0.005. The insert shows the peak near the 
quantum resonance at k = 2ir for k = 11. 



large range of k values by making assumptions about 
the form of Do(n). For example, we can assume that 
D o (0) = D {1) = k 2 /4, and that for n > 2 the diffusion 
rate starts at the initial quantum diffusion rate, D q , and 
decays to zero exponentially with a time constant which 
depends on the quantum break time j^[ . The main prob- 
lem is in determining a form that works well near the 
quantum resonance, where the assumption of exponential 
relaxation in Do(n) breaks down and oscillations occur 
in the diffusion rate, as shown in Fig. [5]. 






FIG. 5: Graph showing diffusion rates as a function of kick 
number for (a) h = 2, (b) K = 6, (c) k = 6.28 and (d) k = 6.4, 
with k — 9 and a = 0.005. Notice the inital quasilinear 
behaviour followed by exponential settling for lower k values 
which contrasts with the oscillatory behaviour for k values 
near the quantum resonance peak. 

The behaviour which we observe in the late time rates 
across the quantum to classical transition for the atom 
optics kicked rotor is similar to that observed by Bhat- 
tacharya et al.nM for the DKR with a continuous posi- 



tion measurement, and leads to similar questions about 
the nature of the quantum to classical transition. The ob- 
servation of such results from a real decoherence process 
is very interesting, and our simulations suggest that these 
results should be readily observable in laboratory experi- 
ments. In fact, hints of unusual behaviour in the momen- 
tum diffusion rates as a function of k have already been 
observed in experiments with cold atoms |l7], [ll| [h| . 

We thank Kurt Jacobs for interesting and stimulating 
discussions and Andrew Doherty for providing the com- 
puter source code of Ref. (§1 , which formed the basis for 
our simulations. This work was supported by a grant 
(UOA016) from the Marsden Fund of the Royal Society 
of New Zealand. 

Note added: After the completion of this work we 
learned of a cold atom KR experiment [M. B. d'Arcy 
et al, Phys. Rev. Lett. 87, 074102 (2001)] demonstrat- 
ing enhanced diffusion resonances of the sort considered 
in this paper. 
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